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Abstract We propose a new formalism of quantum subsystems which allows 
to unify the existing and new methods of reduced description of quantum 
systems. The main mathematical ingredients are completely positive maps 
and correlation functions. In this formalism generalized quantum systems 
can be composed and there is a notion of generalized entanglement. Models 
of fermionic and bosonic systems and also quantum systems described by the 
SU(2) symmetry are studied. 



1 Introduction 

The reduced description of many-particle systems in terms of a relatively 
small number of parameters is a main tool in Statistical Physics. There 
exist several schemes of reductions leading to different mathematical struc- 
tures and involving reduced dynamics as well. The theory of quantum open 
systems within the operator algebra formalism offers perhaps the richest ex- 
ample of such a description. 

Historically, the first general and abstract approach to reduced dynamics 
of quantum systems was introduced by Nakajima [1], Zwanzig [2], and Pri- 
gogine [3]. It is called the projection technique and is based on a projector 
operator V, i.e an operator that satisfies = V. The states of the total 
system are elements of a Banach space B and V projects onto a subspace Bq 
that contains the states of the subsystem, i.e. the reduced system. The pur- 
pose of the projector is to eliminate the irrelevant freedoms of the so-called 
environment, reservoir, or heat bath. This leads for the reduced dynamics 
either to an integro-differential equation with a certain memory kernel or, 
using the time-convolutionless approach [4], to differential equations for the 
states of the subsystem. 
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Another type of reduction can be traced back to the Boltzmann derivation 
of Unear and nonUncar kinetic equations for a gas of particles. Here the 
complete description in terms of A^-particle probability distributions is re- 
placed by a density of particles in the single-particle phase space. In the 
limit of large N and low density one obtains a closed equation of motion. A 
similar approach in quantum mechanics leads to Hartree type equations for 
Hamiltonian dynamics and their extensions to linear and nonlinear quantum 
dynamical semigroups for many-body open systems [5]. 

The aim of this paper is to present a unifying formalism for reduced descrip- 
tions in terms of generalized subsystems (GS). In contrast to the projection 
technique, where only the linear structure of the underlying Banach spaces 
is retained, we heavily use the algebraic structure of quantum theory. The 
projection operator V is replaced by the dual of a completely positive map 
$ from the algebra of observables of the generalized subsystem, often finite 
dimensional, to that of the total system. 

The algebraic formalism allows for a rich structure of the GS's including the 
notion of positivity and the possibility of producing composed GS's with their 
generalized entanglement. This should be compared with a very different 
approach introduced in [6]. The formalism of GS not only unifies several 
known instances of reduced descriptions but it provides also new examples 
like Lie algebraic GS or a quantum like formalism for classical systems. 

The paper is organized as follows: generalized subsystems are introduced 
in Section 2. In Section 3 we show how the GS formalism unifies common 
reduced descriptions such as quantum open systems, coarse-graining, quasi- 
free boson and fermion systems, and mean-field models. Section 4 deals with 
two less common examples: systems that come with a SU(2) symmetry and 
a reduced description of the Koopman formalism [7]. Finally, composition 
and entanglement in space and time are briefly introduced in Section 5. 

2 Generalized subsystems 

We assume in general that the total, usually complex, many-body system is 
described by a unital C*-algebra A its hcrmitian elements corresponding to 
bounded observables. A generalized subsystem is determined by a linearly 
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independent family 

^ {vi,V2,...,Vd) (1) 
of elements of A called a partition. Introducing the element 

d 
i 

we shall distinguish between partitions of unity where m — 1 and general 
ones. 

A partition V generates a completely positive map $ from the algebra of 
complex matrices of dimension d, denoted by A^^ to ^ 

Note that $ is unity preserving if and only if m = 1. The set of states, i.e. 
the linear, positive and normalized functionals on A, is denoted by S{A). We 
consider the pull-back map $* from S{A) to the set of positive functionals 
on Md 

$*(cc;) = o $. (4) 
We can identify the functional $*(a;) with a positive dxd matrix D'^ through 

^*{uj){A)^Tr{D''A), A e Md (5) 

and therefore view as an affine map from S{A) to Md- The matrix 
is called a correlation matrix and its entries are easily computed 

Dr,=ujiv*v,). (6) 

The image of is called a reduced state space and denoted by S{A,\/). A 
reduced state space is easily seen to be a closed convex subset of M'^. In 
particular, the extreme points of S {A, V) are images of pure states a; on ^ 
but generally the converse is not true. The detailed geometrical structure of 
S{A, V) is one of the problems which should be solved for particular examples. 

Two remarks are in order: in some cases the assumption that ^ is a C*- 
algebra can be lifted and partitions in unbounded elements can be considered. 
This leads to domain problems, e.g. one must choose a suitable Hilbert space 
representation of the global system and a subset of states lu for which D'^ 
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exist. Similarly, with some mathematical care, one can extend the definition 
of GS beyond finite partitions, allowing for countably or even continuous 
partitions. 

The map $* is generally not one to one so that infinitely many states on A are 
mapped on a same D, leading to a proper reduction in the description. Still, 
the states in S{A, V) can sometimes encode to a high degree of approximation 
the state of the global system. This can be modelled by a embedding map 

* : S{A, V) ^ S{A) (7) 

which is even possibly non-linear. Consistency is expressed by the require- 
ment 

o * = id . (8) 

A reasonable basis for the choice of an embedding is the maximal entropy 
principle, see [8]: among all states a; on ^ that return a given D, i.e. such 
that $*(cj) = D, we choose for ^(-D) the state of maximal entropy. This 
presupposes both the existence of an entropy S on the global algebra and the 
uniqueness of the constrained maximum. Suppose e.g. that 

S(a;) = -Tra;loga; (9) 

then 

^{D) = Z{a)-' exp{-J2 ^^j^I^j] (10) 

where 

Z(a) = IYexp{- } (11) 

and where a is chosen in such a way that 

= (Z(a)-i exp{- a,jv*vj}) . (12) 

An embedding map ^ allows to construct a reduced dynamics of the GS. 
E.g. if q; = {at 1 1 € R} is the Heisenberg evolution of the global system then 
we can evolve a D{0) e S{A, V) as 

D{t)^^*{-^{D{0))oat). (13) 
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There is no reason to expect a simple dynamical equation for D{t). Neverthe- 
less, simple closed differential equations can be obtained as limiting cases, 
scaling properly the environment, the map ^ and the evolution a. Well- 
known examples are the Markovian semigroup evolution obtained through 
the weak-coupling limit, Hartree-Fock equations, non-linear mean-field equa- 
tions, . . . 

3 Common generalized subsystems 

In this section we briefly rephrase some well-known reduced descriptions of 
quantum systems in terms of GS. 

3.1 Open quantum systems 

The Hilbert space of the total system is the tensor product Tis ® of the 
Hilbert space Tig of the small system of interest and the Hilbert space of the 
environment He- The algebra A of global observables is that of all bounded 
operators of the total system, A = B(Hs ® He)- To a given state to of the 
total system a reduced density matrix D'^ of the subsystem S is assigned 
through the partial trace 

D'^^Tteco. (14) 

Here the state u is identified with its corresponding density matrix which is 
still denoted by the same symbol oo. 

This reduced picture can easily be handled in terms of GS by introducing a 
partition (^1,^2,...) with 

Vj^\<p){j\®lE. (15) 

Here {\j)} is an orthonormal basis in Hs and (f e Hs is an arbitrary nor- 
malized vector. The corresponding map $ is given by 

^A)=Y,Aijv:vj = A(E)lE (16) 

ij 

and $* is the partial trace with respect to the environment 

<^*{uj){A) ^LU {A ^1e), AeB{ns). (17) 
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Obviously, for this example the reduced state space S{A, V) consists of all 
density matrices on Tis. It is also clear that $* is not injective as there are 
many ways to extend a state on <B(7Ys) to the total system. 

In order to obtain a well-defined reduced dynamics one starts by extending 
an arbitrary state D of the system to a global state uj = D ® a;E, i.e. 

*(L>(0)) =L>(0)®u;e. (18) 

Here u}-^ is a suitably chosen reference state of the environment. Obviously 
$* o ^ = id. This embedding of states of the system in global states yields 
a reduced dynamics for the states of the system 

D{Q) ^ D{t) = ^*{U{t)-^{D{0))U{ty) (19) 
^TrE(u{t){D{0)®ujE)U*{t)'^ (20) 

where {[/(t) : i e R} is the reversible evolution of the global system. This 
reduced dynamics is generally very complicated and highly non-Markovian. 
However, in the regime of weak interaction between system and environment 
it is governed by a Markovian master equation of standard form [9] 

±D^-Z[H,D] + ^Y1 ([^" ' ^^a] + [^"^ ' ^a]) • (21) 

a 

This scheme has many applications in various fields of physics including quan- 
tum information processing in the presence of noise. 



3.2 Coarse graining 

Quite often one is only interested in the occupation probabilities of certain 
energy levels or groups of almost degenerate energy levels. These are gener- 
ated by a family of orthogonal projectors (Pi, P2, . . .) with 

P* = P„ PiPj = SijPi, and J2 = ^- (22) 

j 

The corresponding probabilities are 

p,. = Tr(DP,) (23) 
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where D is the reduced density matrix of the system. 

Those probabihties can also be described by a GS defined by a partition 
[Pi ® 1e, P2 <8) 1e, ■ ■ ■)■ The correlation matrices are now always diagonal 
and given by 

D^j = Pi^ij with pj = (jj{Pj ® 1e) . (24) 

Again, under certain assumptions and using a Markovian approximation, one 
can derive the Pauli master equation for the probabilities 

k 

Here aj^ > are transition probabilities per unit time typically computed in 
terms of time-dependent perturbation theory, e.g. Fermi's Golden Rule. 

3.3 One particle description for fermions and bosons 

Assume that the system S is not small but consists of many particles, fermions 
or bosons, described annihilation and creation operators associated to a single 
particle orthonormal basis They satisfy canonical anticommutation or 

commutation relations 

OfcO^ ± a^Qk — 5ki and OfcO^ ± a^Ofc = 0. (26) 

Now the reduction to the single particle description is possible if only additive 
observables of the form 

6 := ^ Bkeala^, B = [B^] Hermitian (27) 

e,k 

are relevant. Therefore, instead of the many-particle density matrix p on 
fermionic or bosonic Fock space J^± the 1-particle density matrix Q — [Qke] 
is used 

Qke = TT{palae) or Tr(g5) = Tr^Jpfe). (28) 

This 1-particle density matrix Q is positive but normalized to the average 
number of particles in the system and not to 1. 
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This reduction can again be phrased in terms of a correlation matrix choosing 
the elements of the partition as 

Vj = aj. (29) 

I.e. 

$(5) = J2Bkialai and $*(a;) = [u;(a*a,)] =: [{i , Qu^j)]- (30) 

The operator is called a symbol, it is positive semidefinite and satisfies 
additionally < 1 for fermions. There is a natural but non-linear map ^ 
from the symbol space to the state space of the full system 

*(Q) = COQ (31) 
where ujq is either the fermionic or the bosonic quasi-free state with symbol 

Q 

(K ■ ■ ■ ■ ■ ■ %i) = de*± ( [Oa , Qib)] ) (32) 

with det_ equal to the permanent. It is again immediate that $*o^' = id. As 
is the general open quantum system setting, the reduced dynamics is quite 
complicated. 

Quite often, however, a system can be well-modelled by essentially noniter- 
acting quasi-particles and the leading dissipative effects are well-approxima- 
ted by processes of quasi-particle decay and production. In these cases the 
Markovian master equation (21) possesses a particularly simple form [9] 

+ 2 ^{ikeiiak , pa*e] + Kp, a}]) + Kke{[al , pae] + [alp, a^Dj 

(33) 

where H = '^^^kdlak with decay matrix 7 = [7^^] > and production ma- 
trix K = [km] > 0. The solution of the master equation (33) is a quasi-free 
dynamical semigroup [9]. This description can e.g. be used to deal with the 
following situations: 

a) decay and production of unstable elementary particles, nuclei, quasi- 
particles, . . . 

b) propagation of a quantum electromagnetic field in media in a linear regime 
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c) transitions between localized electronic states and a large number of low 
lying states accompanied by emission and absorption of energy . . . 

Inserting (33) into (28) one can easily derive a closed evolution equation for 
the 1-particle density matrix Q 

j^Q = -i[h , g] - 1 {(^ ± «:) , g} + «: (34) 

where h := Y^k^k\k){k\. 

Note that for bosons the operators Vj — Uj are unbounded and, moreover, in 
both the bosonic and fermionic case the set of indices {j} may be infinite, 
see the remarks at the end of Section 2. 



3.4 Mean field models 

The simplest setting is that of systems of N identical but distinguishable 
particles. The mean field approximation relies on permutation symmetry 
instead of the usual translation symmetry. Because of this huge symmetry 
group the set of invariant states becomes quite small. In the limit of large N 
the set of permutation invariant states, called exchangeable states, reduces to 
mixtures of permutation invariant product states. Therefore an exchangeable 
state can be seen as a probability measure on the density matrices of a single 
particle. This is called de Finetti's theorem in the classical case and St0rmer's 
theorem for quantum systems. The reduction map is defined by 

$(A) = -^(A01---(8)H \-l<S) ■ ■ ■ <^1<^ A) (35) 

where A is a one-particle observable. It is easily seen that the reduced state 
space is just the full set of single particle states and that the maximal entropy 
embedding is given by 

^(L)) = D®D®---(g)L'. (36) 
The dynamics generated by an A'"-particle Hamiltonian of the form 

i=l i^j 
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will preserve the permutation symmetry of initial states. Here h^^^ and h^"^^ 
are Hcrmitian one and two particle interactions. It then follows that the 
reduced dynamics is described by a non-linear evolution equation in the limit 
of large N 



where Tr2 denotes the partial trace over the second factor in 7i C>?> Ti. being 
the one-particle space. General Markovian dynamics can be handled in a 
similar way. 

4 More examples of generalized subsystems 

We discuss now two classes of GS which go beyond the standard schemes 
of reduced description: Markovian open quantum systems that come with 
a representation of a Lie algebra and a quantum like picture of classical 
systems. 

4.1 Lie algebraic open systems 

Consider a quantum open system containing a Lie algebra Al of operators 
spanned by basis elements Xm = X^ satisfying the commutation relations 



The operators Xj define a partition V = {vj = Xj ® 1e}. 

Assume that the dynamics of the open system is governed by a standard 
Markovian master equation (21) with Lq, = L* . If both the Hamiltonian H 
and the operators belong to Al then the Heisenberg equation of motion 
for a product XmXn reads 




(38) 




(39) 



k 



d_ 

~dt 



{XmXr,) = lXm[H , Xn] + t[Xm , H]X^ " 3 E (t^ 

[La , ^mW^i 



n 



a 




(41) 



(40) 
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This yields a closed evolution equation for the correlation matrix 

^^nm = ^ a{mn; k£)Dek- (42) 

4.2 Angular momentum spaces 

In the following we consider quantum systems with irreducible representa- 
tions of su{2) Lie algebra given by angular momentum operators. We shall 
denote by J^^^ — {Ji\j2\jf^) the irreducible spin £ representation of 
the three standard generators of the rotation group where £ takes values in 
{0,|,1,...}. The jf^ are matrices of dimension 2£ + 1 which satisfy the 
relations 

(^f )* = [J^^ , Jf] = i ^a/4'^ JW • J(^) = l{l + 1)1. (43) 

Here, €^13^ is the totally antisymmetric tensor with 612^ = 1. The Ja^ are 
for a given £ up to a unitary transformation uniquely determined by the 

relations (43) and we shall use the standard convention that jf^ is diagonal 
in the standard basis of C^^+^. To obtain operational partitions of unity we 
have to renormalize the generators 

jW= , ^ - TW (44) 

A reduced state description in terms of jt^^^ can be practically useful for £ ^ 
1, however, the complete description of its structure is here only illustrated 
for the simplest cases £— 1/2, 1, 00. 



4.2.1 Spin 1/2 

The normalized spin 1/2 generators have the form 

and = 





"0 1" 




'0 -i 




1 




i 



1 

7! 



1 
-1 



(45) 



The corresponding completely positive map see (3), sends a matrix A of 
dimension 3 into a matrix of dimension 2 

3 

$(A)= ^ A„^JaJ;3. (46) 
a/3=l 
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The following characterization of positive semi-definite matrices will be use- 
ful. Let A be a square matrix of dimension d and denote by (Ai, A2, . . . , A^^) its 
eigenvalue list repeated according to algebraic multiplicities. The elementary 
symmetric invariant of order k is given by 



= n 

Ac{l,2,...,(i} jeA 
#(A)=fc 



(47) 



Then A is positive semi-definite if and only if ^4 = ^4* and all ek{A) > for 
A; = 1,2, 



The reduced state space consists of all matrices D that satisfy the following 
condition 

TrD^l and Tr{DA) > whenever > 0. (48) 

Note that a D which satisfies (48) is automatically a density matrix, indeed 
$(A) > whenever ^ > because $ is CP. 

One can verify, using the characterization of positivity given in (esi), that 
^{A) e A^2 is positive if and only if ^4 e A^s is of the form 



xi 03 -I- i\i2 a2 + iXi3 
03 + ^^21 X2 ai + iX23 



(49) 



with 



flj G [R, G C and Xi + X2 + X3> 0, 

Xij e R and (xi + X2 + x^f > (Xu - Aai)^ + (A13 - A3i)^ + (A23 
Imposing that Ty[AD) > for all such A implies that 



■ A32)^ 



D 



1 
3 



1 
3 

ia2 —ia\ 



-ia.2 

iai 
1 

3 . 



(50) 



where a := {a-i,a2,as) G IR^ satisfies ||a|| < |. So the reduced state space 
in this example is a ball in and hence affinely isomorphic to the state 
space of a qubit (the Bloch ball). The extreme points are parametrized by 
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||a|| = 1/3, the corresponding density matrices have eigenvalues (2/3, 1/3, 0) 
and are therefore not pure. 



The explicit form of is 



1 —ixs ixi 
ix3 1 iX2 
—ixi ix2 1 



(51) 



where the qubit density matrix p is expanded in Bloch notation 



p = - (1 + cc • cr), cc = (xi,X2,X3) e [R^ and ||cc|| < 1. (52) 



As in this example the dimension of the reduced system is actually larger than 
that of the original there cannot exist a map ^ from the reduced system to 
the full system which is a right inverse of 



4.2.2 Spin 1 



The normalized generators of the spin 1 representation define the operational 
partition of unity (ji, j2, Js) with 



1 

1 1 
1 



J2 = - 



-i 

1 -i 
i 



and ,3 = ^ 



1 

0-1 



The corresponding coarse-graining map is 

kl 

From (47) it follows that ^{A) is positive semi-definite if and only if 

A^A* and A<Tr(A). 



(53) 

(54) 



As any state assigns non-negative values to ^{A) with A < Tr(A) and as 
conversely any functional which takes the value 1 on 1 and is non-negative 
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on ^{A) with A < Tr{A) extends to a state, the reduced state space is also 
characterized by the condition 

A^A*SindA< Tr(A) imphes Tr(DA) > 0. 

It is easily seen that this condition is equivalent to 

D density matrix and D < -. 

To any density matrix D with D < | we may associate a density matrix 

b^l-2D. (55) 

The map D ^ D is affine and one to one from the reduced state space to 
the full state space of M.^. In particular, every extreme point of the reduced 
state space corresponds to a pure state. So, every extreme element of the 
reduced state space is of the form 

^(1-^), (56) 
where P is a one-dimensional projector in AI3. 



4.2.3 Infinite spin 

The limiting operational partition when the total angular momentum tends 
to infinity is given by the relations 

3* = 3h [h : J2] = and cycl. perm., and + jl + j| = 1. (57) 

The Abelian C*-algebra generated by the ji is just the algebra of continuous 
complex- valued functions on the unit sphere in V? . An explicit isomorphism 
is given by 

ji(Q) = sin^cos</7, ^2(^1) = sin^sin(/7, and ^3(^1) = cos^. (58) 

Here, Q = (^, (p) is the usual parametrization of a point on the unit sphere 
in by (co)latitude and longitude: < ^ < tt and < </? < 27r. Again we 
introduce the unity preserving map 

^{A):=Y,Mmn. (59) 

u 
14 



where A is a complex 3x3 matrix. Because of the commutation relations, 
$ is no longer injective. ^{A) is positive if and only if 



n^Y,Akejk{^)ji{Q) 



is a positive function. This is equivalent to impose that 



{x, Ax) > 0, for all x e 



(60) 



or, equivalently that ^4 + ^4^ > 0. Writing 



A^l{A + A-) + l{A-A-^) 



we see that 




) + -TtD{A-A-'). 



As this expression must be non-negative whenever A + A^ > we must 

have that D = D^. Hence the reduced state space now consists of all 3 
dimensional density matrices with real entries. The extreme points are the 
pure states on Ai-i, generated by normalized vectors with real entries. 

We can for this example describe the maximal entropy embedding ^. Let 
us denote hj fl = {9, (p) the usual spherical coordinates on the unit sphere 
S2 in and by the normalized invariant surface measure ^ dcp sin 9d9. 
According to general principles, the maximal entropy embedding ^ is given 



where j = (^1,^25 is) and where A is a real symmetric matrix of dimension 
3. An explicit formula for A in terms of D is not available but the problem 
can be simplified. Let R be an orthogonal transformation of R^, then j 1— > 
Rj can be realized by a change of variables Q Q' that preserves the 
uniform measure. This change of integration variables can therefore be used 
to diagonalize A. A simple inspection shows that in this case the matrix in 
the right hand side of (61) is also diagonal. Therefore, up to an orthogonal 
transformation, we must only solve (61) for diagonal D and A. 



by 




(61) 
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4.3 Quantum description of a classical system 

The general construction can be repeated for the case when ^ is a com- 
mutative algebra and therefore isomorphic to algebra of continuous complex 
functions on a phase space Q. The state space S{Q) of the system is now 
the space of probability measures on Q. Hence we have a one to one corre- 
spondence between a positive, normalized functional cu and the probability 
measure ii{dx) 

ooif)^ j ii{dx)f{x), feC{Q). (62) 

The partition V consists of complex linearly independent functions on Q 

y ^ {vi,V2,...,Vd), VjeC{Q). (63) 

For a probability measure // on Q with associated functional cu we obtain the 
correlation matrix with the elements 



Dij =Lu{v*Vi) = / iJ,{dx)vj{x)vi{x). 
Jn 

Introducing the standard basis {cj} in we can write 

d 

v{x) := '^^Vj{x) Cj, x e f2 

and conclude that 

D"" = I n{dx) \v{x)){v{x)\. (64) 

JQ. 

Therefore the extreme boundary of the reduced state space consists of a 
closed subset of the rank one positive matrices in A^^, pure states in the 
case of a partition of unity. Moreover, any closed convex subset with such a 
boundary can be realized by suitable limits of choices of partitions as in (63). 

5 Composed generalized subsystems 

The formalism for generalized subsystems that was presented here extends 
naturally to composite systems, they are described by higher rank correlation 
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matrices. This has already been used in the context of quantum dynamical 

entropy and quantum symbolic dynamics [10]. If two systems are described 
by the partitions {fQ,}^^;^ and {wk}]^^i, then the composed system is de- 
scribed by the partition {vaWk}aji:- The elements of a correlation matrix of 
the composed system are given by 

Dak;f3£ = (w^V>/3W^) . (65) 

5.1 Generalized entanglement 

When talking about composed systems questions about entanglement natu- 
rally arise. Both the notions of "generahzed entanglement" and "generahzed 
subsystems" have recently emerged in the hterature [6, 11, 12]. The pro- 
posed schemes also deal with projections of states of a large system on a low 
dimensional spaces but the mathematical structures that have been consid- 
ered are not so rich as these presented here. E.g. neither the order structure 
nor the idea of composition is natural in those schemes. 

The basic notion in the conventional approach to entanglement is that of 
separability. We say that the n-party correlation matrix D is separable if it 
can be represented as a classical correlation matrix of the form 

^ai,...,a„,/3i,-,/3n = / ^^idx) fai{x)fp,{x) ■ ■ ■ Un{x)Un{x), (66) 
Jo. 

where is a probability measure on and fa are measurable functions. Using 
the representation (64) it is not difficult to show that the above definition is, 
up to normalization, equivalent to the standard one 

X{ki,k2,...,kn)Pkr^Pk,^---^Pkr,, (67) 

where A(A;i, A;2, . . . , kn) > and where the P^'s are one dimensional projec- 
tors. 

One should notice that for a single system, i.e. for Dag there always exists 
a classical representation (66). Indeed as Q one can take the manifold of all 
normalized vectors ip in the Hilbert and put faiip) = {^a , '4') where {ca} 
is a basis in C^. 
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An entangled bipcirtite correlation matrix 

Consider bosonic systems composed of two parts with corresponding sets of 
annihilation and creation operators (afc,a|.) and {ba,bl^), k,a = 1,2. The 
correlation matrix has the form 

-Dfea/^ = Tr(a;6{ajaa6fe). (68) 

Consider a two-boson state 

m = J2^Mblal\0)®\0), 

u 

with the normalization condition 



where G denotes the 2x2 matrix [7^^] . 
matrix (68) has the form 



and partial transposition of (69) yields 



The eigenvalues of are given by 



For such a state the correlation 



7ii7i2 711721 711722 



'I711I 

712711 I712P 7120^ 7120^ 

72i7n 72i7i2 I721P 721722 

722711 722712 7220^ I722P 



(69) 



7ii| 7i27ii 

7ii7i2 I712P 

721711 722711 

721712 722712 



7iiW 712721 

711722 712722 

I721P 722721 

l2ll^ I722P 



(70) 



±det(|G|) , -(l±Vl-4(det(|G|))2) 

where \G\ — VGG^. So, for det(G) 7^ 0, one eigenvalue of Q^®^ is negative 
and, according to the criterion of partial transposition [13, 14], this means 
that the correlation matrix fl is entangled. 
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5.2 Entanglement in time 



This idea appeared for the first time in the paper by Leggett and Garg [15] 
who proposed to check Beh's inequaUty for correlations corresponding to 
projective measurements at different times. The problem of a "quantumness" 
test for temporal correlations can be easily formulated in our language of 
correlation matrices. If the evolution of the system from time to t is 
described in Heisenberg picture by the completely positive unity preserving 
map At, then we can define a time-dependent correlation function, normalized 
for partitions of unity, by 



One can now apply well-known criteria of separability or measures of entan- 
glement to describe the evolution of "quantumness" encoded in the correla- 
tion matrices of the single system. Examples will be discussed in a future 
publication. One should notice that this approach is related to the formal- 
ism of thermal Green functions in statistical mechanics or quantum field 
theory [16]. In those cases the reference state u; is either thermal or the 
vacuum. 

6 Conclusions 

In this paper we introduced and examined generalized subsystems (GS) as a 
unifying formalism for the reduced description of complex and open quantum 
systems. It covers a number of known examples, like the standard approach 
to open systems with tensor product structure, single particle descriptions 
of many-body systems and Green functions methods. The new examples 
of GS's are quantum systems with symmetries described by Lie algebras. 
Our approach fits well with a large class of approximate evolution equations 
and with the state estimation based on the maximal entropy principle. The 
mathematical formalism involves completely positive maps and correlation 
functions, it has a rich mathematical structure including order relations and 
compositions. This yields a natural notion of generalized entanglement in 
space and time, this issue has only slightly been touched and will be investi- 
gated in the future. 




(71) 
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